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Abstract 

The plasmon oscillations in a cluster of two metallic nanospheres are 
studied theoretically. Particular attention is paid to the case of nearly 
touching spheres. Simple analytical expressions have been found for the 
spectra of plasmon oscillations of different symmetry in this case. A new 
type of the plasmon oscillations, which are strongly localized between the 
spheres, and which totally disappear at separation of the spheres, has 
been discovered. The found plasmon oscillations have a dramatic effect 
on optical properties of an atom localized between the spheres. 

Much attention has been paid, of late, to the experimental and theoreti- 
cal study of optical properties of the metallic nanoparticles. This interest is 
mainly due to a considerable enhancement of local fields near the nanoparti- 
cles. An especially high increase occurs in the case of plasmon polariton reso- 
nances ("plasmons") ^E] or phonon polariton resonances ("phonons")[3]. On 
the basis of this effect one considers quite a number of possible applications. 
One of the most developed is the use of large local fields for enhancement of 
the Raman scattering cross-section 0]. Recent experiments have shown that 
such an increase may achieve 10-14 orders of magnitude, which may help to 
resolve separate molecules [5]- [7]. The local enhancement of the fields can also 
be used to increase the fluorescence intensity and to determine the structure of 
a single DNA strand without using the fluorescent labels [5] |U] ■ By using the 
nanoparticles of complex configuration one can provide enhancement of both 
the absorption and the emission of light by natural and artificial fluorophores 
jl(J|. Of particular interest and promise are the studies of optical properties 
of the clusters of two and more metallic nanoparticles, because by changing 
the cluster's geometry one can effectively control the spectra of the plasmon 
oscillations. This effect makes it possible to produce, for example, new types 
of biosensors ^^"El- A whole series of experimental dl"EI an d theoretical 
|18 | -|27 | studies have been devoted to the two-particle clusters. 
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Among the theoretical works one should note the paper where the 

method for defining the plasmon oscillations of a two-sphere cluster on the 
basis of a plasmon hybridization of separate spheres was proposed. 

Most interesting is the case of nearly touching particles where the enhance- 
ment of local fields is particularly high. This is the most complicated case as 
well. In this interesting case was not actually considered, and the paper 

[TO] being especially devoted to that problem, turned out to be misleading. Bas- 
ing on the numerical consideration of one of the particular cases, the conclusion 
was made in j^U] about failure of the Mie theory for the closely set nanosphcres. 

In the present work we will investigate the most interesting and the most 
complicated case of the plasmon (phonon) oscillations in the cluster of two 
almost touching metallic nanospheres. We shall derive analytical expressions 
for the plasmon spectra in this region. It will be shown that in this system 
there exists a new type of plasmon oscillations (plasmonic molecule), which is 
strongly localized between the spheres, and which cannot be found within the 
framework of the hybridization method JS| • We shall also show that new modes 
can be excited only by a localized light source (an excited atom, molecule, or a 
quantum dot) that should lead to a dramatic changing of optical properties of 
an atom. A geometry of the problem is illustrated in Fig.l. In the case of E\ < 0, 
£2 = 1 the geometry corresponds to two metallic nanoparticles, and in the case 
of E\ — 1, £2 < to two spherical cavities (bubbles) in a metal. Hereinafter 
we assume that the dimensions of the studied particles are sufficiently small 
in order one could use the longwavelength perturbation theory [2H| but at the 
same time, they are essentially large for one could neglect the non-local effects 
(space dispersion of permittivity) 2!) . For simplicity here we will consider only 
the case of equal nanospheres. 

To determine the spectrum of the plasmon oscillations in the first order of 
longwavelength perturbation theory |28| one must solve the equations: 

divD = 0, rotE = (1) 

with usual boundary conditions of continuity of the tangential component of E 
and the normal component of D. To find further small corrections to spectrum 
of plasmon oscillations one must use procedure elaborated in |28j . By means of 
a conventional substitution, E = — V$, the system of equations Q is reduced 
to the Laplace equation for the potential which can be conveniently solved in 
a bispherical system of coordinates. The bispherical coordinates (—00 < 77 < 00, 
0<£<7r, < (f < 2ir) are connected with the Cartesian coordinates by the 
relations [3U| 



sin £ cos (f 
cosh T] — cos £ ' 

sin £ sin cp 

cosh T] — cos £ ' 

sinh?] 
t 

cosh 77 — cos£ 



(2) 
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The surface rj = —rjo < is the sphere of the radius Rq = a/ sinhyyo. 
The other sphere can be set in the similar way by the equality 77 = 770 > 0. 
The dimensional constant a is the half distance between the poles of a bi- 
spherical system of coordinates, and is determined by the next expression: 
a = ^/i?, 2 2 /4 — Rq, where R12 is the distance between the centers of the first 
and second spheres. 

The potential in space between the spheres ( — 770 < rj < r/o) can be written 
in the form 



$ = -v/coshry-cosC^ X F ™ ( C0S ^ 

n— m— 

( 0m „c„s M + /3m „sin(„,.,) ) ™h((, 1 + I), ) 

+ {irrm cos (m(p) + S mn sin (m(p)) sinh f f n + 77 j | (3) 

where a mn , (3 m m Jmn, <5mn are the coefficients to be determined. 

For the potential inside the nanospheres we have, respectively, {q < ~t]q < 0) 



oo n 

= -Vcosh?? _ cose X e( n+ ^)"P™ (cosO 

n— m— 

{a™™ co s (rrup) + CI sin } ( 4 ) 



and (ry > 770 > 0) 



$W = -^coshry-cos^ X e _ ( n+ 5)"P™(cosO 

n— m— 

ja£^ cos (my) + b$ n sin (my) } (5) 

in Q and J5j a^L &ron an d a™, are the coefficients to be determined. 

Expressions for the potentials in the form of ©-© have the most general 
form, and make it possible to determine all the modes irrespective of their 
symmetry. 

By using usual boundary conditions at the surface of each of the nanospheres 
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(6) 



and 
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*(»7 = »7o) = (V = Vo) , 
d d 

(v = Vo) = £ i^ $(2) (v = Vo) (7) 

and the recurrent relations for the Legendre functions one can derive an infinite 
system of equations for all the coefficients: a mn , 7 mn , j3 mn , S mn , a mn , and b mn . 
In the most interesting case of the two identical nanospheres these equations 
can be written down in the following form 



t {(n - m) 5m, n -i + (sinh 770 - (2n + 1) cosh 770) a mn 
+ (n + m+l) a m ,n+i} 

-(n-m) tanh - 770^ a m ,n-i 

+ cosh 770 ^tanhr^o + (2n + l)tanh ^fn + «mn 

- (n + m+ 1) tanh [ ( n + | ) % J a m ,n+i (8) 



and 



r {(n - m) 7 m , n -i + (sinh 770 - (2n + 1) cosh 770) j mr , 
+ (n + m+ l)7m,n+l} 

-(n-m) coth ( f n - ^ r/o^ 7 m , n -i 

+ cosh 7/0 ^tanh 770 + (2n + 1) coth ( f n + 770^ ^ 



- (n + m + 1) coth I I n + - I 770 J 7,„,„+i (9) 

where r = ei/e 2 and a mn = a mn cosh ((n + |) 770), 7 m „ = 7 m „ sinh ((n + \) 770) 
Note that the recurrent equations for the second series of the coefficients j3 mn 
and 5 mn have the form JHJl and (the substitutions a mn — > /3 m „ and 7 mn — > 
<5 mn are to be made). 

The coefficients (£1 and amn for the fields inside the nanospheres can be 
found from the relations 



a (i) g-^+l)^ = a 

(2) p -(n+i)no _ ~ , ~ /iq\ 
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7mn * ^mn- 



The infinite systems of homogeneous equations (jSJ) and © will have nonzero 
solution only for special (resonant) values of r = ei/e^- To find these resonant 
values one should solve the generalized eigenvalue problem. Note that the co- 
efficients a mn and 7 mn (which correspond to solution of different parity in rj) 
enter into JHJ and independently. 

Figure 2 illustrates the dependence of the resonant permittivity on distance 
between spheres for m = 1. For any other m we have similar spectra. One can 
see in Fig. 2 that there are three types of the modes: L-, M- and T-types. As 
seen in Fig. 2, if the distance between the nanospheres increases the resonant 
permittivity for the T- and L-modes takes the known [23 E] values r = — H±l ) 
which correspond to the plasmon oscillations of different electric multipolarity 
excited in a single sphere: n = 1 corresponds to the dipole oscillations, n = 2 
corresponds to the quadrupole oscillations and so on. Such modes have been 
well studied and are described by hybridization of modes of the separate spheres 
|18| . As for the M-modes, it seems that their appearance had not been observed 
earlier. These modes have a strongly localized character, because they can exist 
in the clusters with < 1.2 only. Such modes can be considered as bound 
states of plasmons of the different spheres (plasmonic molecule). In order to 
excite such modes one needs a localized source of electromagnetic field. In the 
papers |^] |2Z] studying the surface modes of a cluster one did not observe the 
appearance of such modes because they considered an excitation of a cluster by 
a plane electromagnetic wave which can not excite new M-modes. 

From the above Figure 2 it is seen that the behavior of the resonant permit- 
tivity in the region of the nearly touching spheres becomes very complicated. To 
describe that complex behavior we consider JHJ) at the limit of ?/ — > (|^ — > 1). 
In this case we have found an asymptotic expression for the resonance parameter 
r (relative permittivity). 

For the M-modes we have 



where m = 0, 1, 2, . . .; L = 1, 2, 3, . . . Of course expressions (lllf) and (|12fl are 
in agrement with Fig. 2 in the region — > 1. 

From these expressions it follows that the resonant relative permittivity 
of the cluster r tends to zero (M-modes) or to the infinity (L-modes) if the 
nanospheres are drawing together. Tending of the relative permittivity to zero 
corresponds to two physically different cases, namely, two metallic nanospheres 



T m .M = — = - (M + S m ) arccosh— ^ + . . . 



(11) 




where m = 0, 1, 2, . . .; M = 1, 2, 3, . . .; S = 1/2 and 8 X « 0.914. 
For the L-modes one can obtain the following analytical solution 




(12) 
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(s\ = 0, e-2 = 1) or two cavities in a metal (ei = 1, e% — > — oo). The analogous 
situation occurs when the relative permittivity tends to the infinity. 

The eigenvectors of the L- and M-modes correspond to the resonant per- 
mittivity l|12fl and at the limit of the closely set spheres can be calculated 
analytically as well, 



/ "I 



M 



\ 



M 



\ 



M 






7m,L 



/ 




(13) 



Note that the analogous expression for 5 m: L coincides with expression l|13fl 
for 7m, l, and /3o,m = 0. One can substitute these vectors into ©-© to find 
the expressions for the potential eigenfunctions. 

For the eigenfunctions of the M-mode (m = 0) in space between the spheres 
we obtain 



M-l 



$0,M 



-r^j- \J cosh rj - cos 1 ^2 cosh ( ( n + ~ J rjj P n (cos£) 



+ - ^/cosh T] - cos | cosh ^ I M +-)»?) -Pm (cos f ) ( I 4 ) 

Inside the first sphere we have, correspondingly, 



(i) 

0,M 



170^0 



M-l 



£ (l + (»+ 5) ^oj e( n+ i>'P n (cost) 

+ iv/coshr/-cose ^1 + fjlf + 770^) ei M+ ^ r 'P M (cos£) (15) 

Inside the second sphere we obtain the analogous expression. 

For the eigenfunctions of the i^mode in space between the spheres and inside 
the first sphere we have, correspondingly, 



m,L\, lo _ 



— y/2 (cosh rj — cos£) sinh ^ (^L — 77 
Pl-i (cos£) sin fm<£ + 



(16) 



G 



and 



7) ->o a V 2 



pm 



^(cosh^-cosOe^ - *)" 
x (cos sin (m</?+ (17) 



Inside the second sphere we get the analogous expression. 

Knowing the eigenfunctions one can calculate a surface charge a^' that is 
stored on each of the nanospheres (i=l, 2) 



,., 7- — 1 . <9$ w 

ct w = — (cosh 770 - cos£) (-1)' 



4na drj 



(18) 



V = Vi 



Using the asymptotic expressions for the potentials l|14 |l -l|17 jl one can write 
down the expression for the charge for the M-mode (m = 0) in the form 



T (i) _(a) ^ (1-cosQ 



3/2 



1 M_1 \ 

- Y, ( 2n + 1 ) p n (cos + (2M + 1) P M (cos (19) 

n=0 / 



\ M 

For the i^mode we get, respectively, 

(1) _ _ (2) _ (l-cosQ 3/2 (2L-1) 2 
r m,L ° m ,L~ 4^/2^2 2(m + L)- 



T Pr_ 1 (cose)sin(m^ + J) (20) 



We did not succeed in finding the analytical solution for the T-mode. More- 
over, it is difficult to find limit at — * 1 of these curves. To investigate this 
limit let us consider the limiting case of the touching spheres = 1) within 
the coordinate system of the touching spheres {31 ] . In this case, the spectrum 
of the plasmon oscillations that are symmetrical to the z—0 plane is determined 
by the integral equation 



A' oo 

K m (A') J dXA mX I m (A) + I m (A') J d\A mX K m (A) 

A' 

' (r + tanh(A'))A mA , (21) 



T - 1 

where I m and K m are the modified Bcssel function and the Macdonald function, 
and A m \ is the Hankel transform of the potential eigenfunction. 

The integral equation (|21(l for the eigenvalues r mi T can be solved numerically 
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to,i » -1.696, 7D, a w -1.355, r , 3 « -1.237, 

n,i fa -1.799, Ti,2 « -1.386, r 1)3 » -1.182 (22) 

Thus, in the case of the touching nanospheres the intrinsic resonant values 
of T-modes lie within the interval [-2,-1]. The numerical solution of (|22|l agrees 
well with the functions represented in Fig. 2 in the region 

of i& ~" 1 (T-rnodes). 

In the case of antisymmetrical (longitudinal) Z^modes the integral equation 
has the form: 

A' oo 

K m (A') j d\A mX I m (A) + I m (A') J d\A mX K m (A) 

A' 

= — ^— (r + coth (A')) A m \i (23) 

T — 1 

which seems to have only trivial solutions at finite Ex/e2- This result agrees with 
the asymptotic of (|12|) and disagrees with a numerical solution found in |2(Jj . 
Perhaps, its difference from [20] is due to insufficient accuracy of the calculations 
in H3. 

Figure 3 illustrates surface charge distributions derived from an accurate 
numerical solution of the recurrent equations (JSJ and © for m = 1. For any 
other m we have the similar distributions. To plot these distributions we have 
used a polar system of coordinates (p, 9) for each sphere in plane y — 0. The 
surface charge distribution was depicted with p = Ro (1 + a iff)) polar curve 
where a (9) is the dimensionlcss surface charge density as a function of polar 
angle of corresponding sphere. The M-mode and its vividly local character are 
of most interest. This mode is actually localized in the region with characteristic 
size of the order of size of the gap between the nanospheres. The T- and L- 
modes are due to the charge oscillations over the whole surface of spheres and 
are analogous to the charge distribution on separate spheres. Note that the 
L-modes have also appreciable localization between the spheres. 

It is very important that new M- modes can not be excited by uniform 
electric field. To excite such modes one must use nonuniform field such as 
dipole radiation of an atom or molecule placed near the cluster. On the other 
hand, the excited plasmon modes will influence on decay rate and fluorescence 
of an atom or molecule. The modification of decay rate by new modes can be 
used to detect them. 

The radiative decay rate ry radlatme f an atom placed near any nanobody 
can be described by j33 ESI 



-.radiative I J |2 
7 \0-total\ 



(24) 



7o do- 
where 70 is the atom spontaneous decay rate in free space; do, d tota ; = do + (5d , 
the dipole transition moment and the total dipole moment of an atom + cluster, 
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respectively. In our case the expression for the dipole moment induced in the 
cluster takes the form pH] 



oo 



Sd 



•X 



n (n + 1) (doV ) 



oo 



Sd 



■y 



%/2a^n(n+l)(d V')/3: 



71=1 



OO 



Sd z 



V2a ( 2n + 1) (do V) 7a 



(25) 



where V' means gradient over the atom's coordinates; a\ n , (3\ n , 7o n are the 
coefficients in found with taking into account the atomic dipole source. 

Figure 4 illustrates the dependence of the spontaneous decay rate of an atom 
near two nanospheres (|24|) . I|25|) on a transition wavelength. The spheres are 
placed in vacuum (£2 = 1) and made of SiC, where phonon resonances can be 
excited in the infrared [35]. Strictly speaking SiC is uniaxial hexagonal crystal 
with anisotropic tensor of dielectric permittivity. However, strong TO_L (797 
cm -1 ) and TO|| (788 cm' 1 ) modes have very similar dielectric properties. It 
allows us to consider SiC crystal in the first approximation as isotropic one in 
the region of 788-797 cm -1 . In Fig. 4 the influence of these modes on sponta- 
neous emission of an atom is shown. As seen in Fig. 4(a), the T- and M-modes 
can be excited by the atom with the transversal dipole (the atom lies precisely 
between the spheres). The excitation of only these modes can be explained by 
the symmetry of an induced charge at the surface of the spheres, in respect to 
the plane with z = 0. As the real part of the permittivity increases (decrease 
in the transition wavelength) then T-modes are replaced by M-modes, in ac- 
cordance with the functions of Fig. 2, at transition over the E\ = —1 point. It 
is also important to mention that the spontaneous decay rate for the M-modes 
is substantially higher than that for the T-modes. In Fig. 4(b) one can observe 
only the Z^modes in the spontaneous decay rate spectrum, provided that the 
dipole moment of an atom is oriented along z-axis (the atom lies precisely be- 
tween the spheres) . This fact is due to the symmetry of the excitation source 
also, i.e. the charge distribution on spheres is to be antisymmetrical to the plane 
with z = 0. In the same Figure one can clearly see an increase in the number 
of the resonances as real part of the permittivity approaches E\ = — 1, which 
agrees with the functions represented in Fig. 2. As the transition wavelength 
decreases (going over e± = —1) no resonant phenomena are observed for such 
an orientation of the dipole, except for the dip at e± = 0. 

Thus, in this paper the spectrum of plasmon oscillations occurring in a 
cluster of two closely set nanospheres, or spherical bubbles, has been considered 
in detail, and the respective analytical expressions have been found. A new type 
of strongly localized plasmons oscillations has been discovered. Such modes can 
be considered as bound states of plasmons of the different spheres (plasmonic 
molecule) . New modes exert influence on the optical properties of an atom near 
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a gap between spheres, and can be used to control the decay rates of an atom 
and molecule. From the other hand well pronounced dependencies of an atom 
optical properties on distance between spheres, dipole moment orientation and 
permittivity can be used for elaboration of new types of optical nanosensors. 
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List of Figure Captions 



Fig.l Geometry of the problem. 

Fig. 2 Resonant permittivity as the function of normalized distance between 
two identical spheres (m = 1). 

Fig. 3 Surface charge density (in the relative units) for different modes with 
m = 1 (the plane y = is shown), (a) L- modes; (b) T-modes; (c) M-modes. 
The distance between spheres is = 1.1. The broken line denotes a surface 
of the sphere; the colored lines correspond to the mode indices: red-1; green-2; 
bluc-3 (cf. Fig.2). 

Fig. 4 The spontaneous emission radiative decay rate of an atom near the 
two nanosphcres made of SiC. The spheres have the radii of 50 nm; the distance 
between their centers is of 101 nm. The atom is placed precisely between the 
spheres and its dipole moment is directed: (a) transversally to the rotation axis; 
(b) longitudinally to the rotation axis. 
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Figure 2: 
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